arXiv:1501.06637vl [quant-ph] 27 Jan 2015 


Algebraic structure of the two-qubit quantum Rabi model 
and its solvability using Bogoliubov operators 


Jie Peng 12 , Zhongzhou Ren 2 3 4 5 , Haitao Yang 6 , Guangjie Guo 7 , Xin 
Zhang 2 3 , Guoxing Ju 2 , Xiaoyong Guo 8 , Chaosheng Deng 1 , Guolin Hao 1 

laboratory for Quantum Engineering and Micro-Nano Energy Technology and School of 
Physics and Optoelectronics, Xiangtan University, Hunan 411105, China 
2 Key Laboratory of Modern Acoustics and Department of Physics, Nanjing University, 
Nanjing 210093, China 

3 Joint Center of Nuclear Science and Technology, Nanjing University, Nanjing 210093, 

China 

4 Center of Theoretical Nuclear Physics, National Laboratory of Heavy-Ion Accelerator, 
Lanzhou 730000, China 

5 Kavli Institute for Theoretical Physics China, Beijing 100190, China 

6 College of physics and electronic information engineering, Zhaotong University, Zhaotong 

657000, China 

1 Department of Physics, Xingtai University, Xingtai 054001, China 

8 School of Science, Tianjin University of Science and Technology, Tianjin 300457, China 

E-mail: pengjiel45@163 . com, zren@nju.edu.cn and jugx@nju.edu.cn 

Abstract. We have found the algebraic structure of the two-qubit quantum Rabi model 
behind the possibility of its novel quasi-exact solutions with finite photon numbers by 
analyzing the Hamiltonian in the photon number space. The quasi-exact eigenstates with at 
most 1 photon exist in the whole qubit-photon coupling regime with constant eigenenergy 
equal to single photon energy Hui, which can be clear demonstrated from the Hamiltonian 
structure. With similar method, we find these special “dark states"-1 ike eigenstates commonly 
exist for the two-qubit Jaynes-Cummings model, with E = NHuj (N = —1,0,1,...), and one 
of them is also the eigenstate of the two-qubit quantum Rabi model, which may provide some 
interesting application in a simper way. Besides, using Bogoliubov operators, we analytically 
retrieve the solution of the general two-qubit quantum Rabi model. In this more concise 
and physical way, without using Bargmann space, we clearly see how the eigenvalues of 
the infinite-dimensional two-qubit quantum Rabi Hamiltonian are determined by convergent 
power series, so that the solution can reach arbitrary accuracy reasonably because of the 
convergence property. 
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1. Introduction 

The quantum Rabi model fU describes the interaction between a bosonic mode and a two- 
level system—probably the simplest interaction between light and matter. Its semiclassical 
form was first introduced by Rabi in nuclear magnetic resonance 0. In 1963, Jaynes 
and Cummings 01 found its application in describing the interaction between a two-level 
molecular and a single mode photon field. With the developments of experiments, many 
systems can be described by this model in quantum optics 0, condensed matter 0], cavity 
quantum electrodynamics (QED) Q, circuit QED 10, quantum dots 01, trapped ions [(8j and 
so on. Although this model takes a very simple form, its analytical solution was not so easy 
to obtain, so various approximations were made, one of which is the famous “rotating-wave 
approximation” 0. In 2011, its solution was analytically found by Braak [J9]] in the Bargmann 
space lITOil . It can describe the ultrastrong qubit-photon coupling regime, which has been 
reached in recent circuit QED experiments |[TT|. where the “rotating wave approximation” 
breaks down. After that, various researches are done to the full Rabi Hamiltonian, including 
recovering the solution of the Rabi model Ifl2] - [l4ll . real-time dynamics lfl5l . the solution of the 
two-qubit Rabi Hamiltonian lTT6l420ll . dynamical correlation functions lITTIl . and so on lf22l - f3Til . 

Two-qubit system is basic and fundamental to the construction of the universal quantum 
gate. Various qubit-qubit interactions are applied to generate qubit-qubit entanglement and 
realize quantum computation 1132113311 . one of which is mediated by a resonant cavity, described 
by the two-qubit quantum Rabi model llT9ll . In this case, the ultrafast two-qubit quantum gate 
can be constructed in the ultrafast qubit-photon coupling regime ll34l . Besides, the distant 
qubits can be coupled through a resonant cavity and the coherent quantum state storage 
and transfer can be realized f[35ll . Working for the whole qubit-photon coupling regime, the 
two-qubit quantum Rabi model can be applied in many systems in quantum optics lf36l and 
quantum information lf37ll . Its analytical solution was obtained in llT9l by means of Bargmann 
space approach, and also in ll20ll with extended coherent states representation. One interesting 
result is that there exist coupling-dependent eigenstates in the whole coupling regime with 
constant eigenenergy-reminiscent of “dark states”, but they are coupling-dependent and the 
photon number is bounded from above at 1, which is novel and interesting. Besides, there are 
quasi-exact solutions with finite photon numbers N, which are not presented in the one-qubit 
Rabi model. These special solutions may have some interesting application, however, the 
algebraic structure behind the possibility of these special solutions needs to be clarified. 

In this paper, we have clarified the algebraic structure of the two-qubit quantum Rabi 
model for its special quasi-exact solutions with finite photon numbers found in lfl9l . By 
analyzing the Hamiltonian in the photon number space, we find the condition for closed 
subspace, i.e. the algebraic structure are related with the permutation symmetry of the qubit- 
photon coupling terms for the two qubits. Even more interestingly, the quasi-exact solution 
with at most 1 photon exists in the whole coupling regime with constant eigenenergy equal 
to single photon energy hu, which can be clearly found from the algebraic structure. These 
eigenstates are partly like “dark states”, but are coupling dependent and the photon number is 
bounded from above, so they may have some interesting. According to the algebraic structure 
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of the two-qubit quantum Rabi model, we may conjecture there are similar “dark states”- 
like solutions to those models with homogenous qubits-photon coupling terms. For example, 
we consider the two-qubit Jaynes-Cummings model If35l . which is commonly applied for 
simplicity in the weak coupling regime f[38ll . Very interestingly, under similar condition, we 
find many “dark states”-like eigenstates, existing in the whole coupling regime with constant 
eigenenergy E = Nhu (N = —1,0,1,...), one of which is also the eigenstate of the 
two-qubit Rabi model. Since the Jaynes-Cumming model is simper than the Rabi model, 
these eigenstates may provide some interesting application easier. On the other hand, we 
analytically retrieve the solution of the two-qubit quantum Rabi model, using Bogoliubov 
operators. With this more physical and straightforward method, we find a way to obtain 
its solution by convergent power series, so that we can make reasonable cutoff in practical 
calculation and the solution can reach arbitrary accuracy. 

The paper is organized as follows. In section El we clarify the algebraic structure behind 
the possibility of quasi-exact solutions with finite photon numbers obtained in lfl9ll and also 
find the special “dark states”-like solutions of the two-qubit Jaynes-Cummings model. In 
section [3j we analytically retrieve the solution of the two-qubit quantum Rabi model using 
Bogoliubov operators. Finally, we make some conclusions in section|4l 


2. Algebraic structure for quasi-exact solutions with finite photon numbers 


The Hamiltonian of the two-qubit quantum Rabi model reads (h = 1) Ifl7lfl9 1 


H tq = l oa ] a + gicr lx (a + a t ) + g 2 cr2x(a + a t ) + Aicr lz + A 2 a 2z , 


( 1 ) 


where and a are the single mode photon creation and annihilation operators with frequency 
u, respectively, a, (i = x, y, z) are the Pauli matrices. 2A l5 2A 2 are the energy level 
splittings of the two qubits. g\ and g 2 are the qubit-photon coupling constants for the two 
qubits respectively. There are quasi-exact solutions with finite photon numbers N obtained by 
analyzing the recurrence relation of the coefficients in lfl9l . However, the algebraic structure 
behind the possibility of these novel exceptional solutions needs to be clarified. 

Quasi-exact solutions with finite photon numbers N correspond to the existence of 
closed subspace in the photon number representation, i.e. the algebraic structure. Here we 
demonstrate the closed subspace are related with the permutation symmetry of the qubit- 
photon coupling terms by analyzing the structure of the Hamiltonian in the photon number 
space. H tq ([[]) process a Z 2 symmetry with the transformation R = exp(f7ra t a) 0 cr lz 0 a 2z . 
Taking odd parity for example, supposing the initial state | if)) is in a subspace formed by 
{| M,e,g), | M,g,e), \M + l,g,g), \M+l,e,e), ■ ■ ■ , \N - l,g,g), \N - l,e,e), \N,e,g), 

| A, g, e)}, with the coefficient {ci jM , c 2 , M , c 1 >M +i, c 2 ,m+i, • • •, c 1>N , c 2)N }, where M and N 
are even, then the Hamiltonian reads (u is set to 1) 


( 0 0 \/Mgi \[Mg 2 0 

0 0 y/Mg 2 \[Mg\ 0 

y/Mgi \[Mg 2 M + A x - A 2 0 y/M + lg\ 

V Mg 2 \JMgi 0 M + A 2 — A x \JM 4- 1 g 2 

\ . 


0 

0 

y/M + lg 2 

VM + lgi 


\ 

/ 
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V 


y/N g\ y/N g 2 N + A\ — A 2 0 

y/Ng2 VN 91 0 N + A 2 — Ai 

0 0 y/N + 1 g \ y/N + lg 2 

0 0 y/N - f- lg 2 y/NTlgi 


... ... \ 

VA + I51 y/N + lg2 
y/N + lg 2 y/N +lgi 
0 0 

0 0 / 


( 2 ) 


If for | ijj') = H | -0), the coefficients of {| A + 1, g, g), |A + l,e,e)} and {\M — 1, g, g), \M— 
1, e, e)} equal to 0, then this subspace is closed. For the first case, we obtain 


VA + lgic^ N + y/N + lg2C2,N — 0, (3) 

VN + lg 2 ci jN + y/N + lgic 2 ,N = 0, (4) 


where c^ N and c 2 ,n are the coefficients of | N, e, g) and \N, g. e ) respectively. From equations 

© and © and < 71 , g 2 > 0, we obtain g\ = g 2 and c^n = ~C 2 ,n- By using the time- 

independent Schodinger equation, we obtain 

V^V^ic^jv-i + y/~N g 2 C2 } N—i + (N + Ai — A 2 )ci j at = (5) 

'/~Ng2C\,N-i + V^.gi c 2 ,n-\ + (N + A 2 — A 1 )c2 ] at = Ec2,n, (6) 


so that 


E = N, (7) 

(A 2 — Ai)ci ; Ar = (VN (/iC^jv-I + V / iV(?2C2,Ar-l). (8) 

For the special case Ai = A 2 and c \ v -1 = c 2 ) n-i = 0, there is a invariant subspace formed 
by {| A, e, g), \N, g , e)}, and the eigenstate is 

|^)tv = ^d-^’ 9, e) - |iV, e, g)), (9) 

which is the famous “dark state” Ill9[l39l . where the spin singlet is decoupled from the photon 
field. 

If Ai ^ A 2 , considering the coefficient of {| M — 1 ,g,g), \M — l,e, e)} must be 
0, it is required that E = M, which contradicts with E = N, so that the only possible 
choice is M = 0. Now we have obtained a closed subspace (algebraic structure) formed by 

{|0, e, g), 10, g, e), |1 ,g,g),\l,e,e), ■ ■ ■ ,\N - 1 ,g,g), \N - l,e,e), | N,e,g), \N,g,e)}, with 
the condition 


9i = 92, (10) 

E = N. (11) 

Then by using the time-independent Schodinger equation, we can obtain quasi-exact solutions 
with finite photon number N for certain choice of parameters Ai, A 2 , and g = gi + f) 2 - For 
example, if A = 2, the determinant of the matrix 


/ Ai —A 2 —2 

0 

5/2 

5/2 

0 

0 

0 

A 2 - Ai - 2 

5/2 

5/2 

0 

0 

g /2 

5/2 

—Ai — A 2 — 1 

0 

y/2g/2 

y/2g/2 

5/2 

5/2 

0 

A 2 + Ai — 1 

V2g/2 

y/2g/2 

0 

0 

y/2g/2 

V2g/2 

Ai — A 2 

0 

V 0 

0 

y/2g/2 

y/2g/2 

0 

A 2 — A 


( 12 ) 
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must equal to 0, which gives 


(A? - A=)[(4 - (A, - A 2 ) 2 )(l - (A, + A 2 ) 2 ) - 2g 2 ] 


0 . 


(13) 


This is the condition for an odd parity solution with photon number bounded from above 
at N = 2, coinciding with lfl9ll . which depends on A 1? A 2 and g. So now, we have found the 
algebraic structure and quasi-exact solutions with finite photon numbers N. Furthermore, it 
is very interesting for the solution with N — 1, whose existing condition is independent of g. 
The closed subspace is formed by {|0, e, g ), |0, g , e), 11, g , g), 11, e, e)}, and the condition is 


Ai — A 2 — 1 

0 

g/ 2 

g/ 2 

0 

A 2 — Ai — 1 

g/ 2 

<7/2 

g/ 2 

g/ 2 

—Ai — A 2 

0 

9/ 2 

g/ 2 

0 

A 2 +A 


which gives 

(A 1 + A 2 ) 2 [(A 1 -A 2 ) 2 -1] = 0, 


(14) 


(15) 


which is independent of g, coinciding with [fl9ll . So for Ai — A 2 = 1 = fku and 
Ai — A 2 = — 1 = —hui, we obtain two quasi-exact solutions 

1^) g i = ^ 2 ( Al + A ^ |o+ |1 ,g,g) - |l,e,e)^ , (16) 

\^)g 2 = ^ 2 ( Al + A ^ |0 ; g y e) + |1 ,g,g) - |l,e,e)^ , (17) 

respectively, where J\f = ^/4(Ai + A 2 ) 2 + 2 g 2 /g. For example, choosing A x = 1.4, 


A 2 = 0.4, gi = g 2 , the numerical spectrum of the two-qubit quantum Rabi model is shown 
in figure [Q The horizontal line at E = 1 = huj corresponds to the special eigenstate 
\ip) g i (equation (fl6l)). This eigenstate exists in the whole coupling regime with constant 
eigenenergy, like “dark states”, but are coupling dependent, and with at most 1 photon. 

For even parity, similarly, we obtain one such special eigenstate 

\^)e = ^ Al ~ A ^ |0,e,e) - 1 1, e, g) + |1 ,g,e) \ , (18) 

with the condition Ai + A 2 = 1 = ui, gi = g 2 and E = 1 = ftw, consistent with [fl9l . Now, 
we have demonstrated all the exceptional eigenstates of the two-qubit quantum Rabi model 
with finite photon numbers presented in lfT9] by finding its algebraic structure in the photon 
number space. 

The special eigenstates \ip) g i, \ip) g 2 an£ l \i>)e originate from the permutation symmetry of 
the qubit-photon coupling terms, and we may conjecture there are similar solutions for similar 
models. In the weak-coupling regime, the Rabi model can reduce to Jaynes-Cummings model 
by the rotating-wave approximation, so if there are similar special eigenstates for the two- 
qubit Jaynes-Cunmmings model, we may find its application in a simpler way. Now we try to 
find similar structure for the two-qubit Jaynes-Cunmmings model |[38l 

H tjc = a) a + gi(af a + cr 1 ~a t ) + g 2 {^ a + ovTa t ) + A x a lz + A 2 cr 2 .. (19) 
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Figure 1 . The numerical spectrum of two-qubit quantum Rabi model with Ai = 1.4, A 2 = 
0.4, u) = 1, gi = g 2, 0 < g = g± + g 2 < 3. E + and E_ are solutions with even and odd 
parity respectively. 


It is easy to find C = a^a + |(oTz + cr 2z + 2) commutes with H t]C , so there is a conserved 
quantity C. Interestingly, \ip) e (equation (ITSl) ) has a conserved quantity C — 2, and it is 
easy to testify \ip) e is also an eigenstate of H t j C existing in the whole coupling regime with 
constant eigenenergy for rq = g 2 and Ai + A 2 = 1. To find out all such kinds of eigenstates, 
we study the eigenproblem of H t j C . For C = N (N > 1 ), the Hamiltonian in the subspace 
{|iV — 2, e, e), | N - l,e,g), \N - l,g,e), \N, g, g)} reads 

( N — 2 + Ai + A 2 yjN — 1 <72 y/N — lg\ 0 ^ 

VN - lg 2 N — 1 + Ai — A 2 0 y/Ngi 

y/N - lgi 0 N - 1 - Ai + A 2 VNg 2 ' { 

\ 0 VNgi VNg 2 N - Ai - A 2 / 

Using the time-independent Schodinger equation, we find the eigenvalues E is deteimined by 

(E — N + Ar + A 2 )(E - N + 2 - A 1 - A 2 )[(E - N + l) 2 - (A 1 - A 2 ) 2 ] 

+ G g\ + g 2 )[(E -N + 1 )(E - N + A 1 + A 2 ) - 2 N(E -N + l) 2 ] 

+ (ai - gi)[(gl - g!)(N 2 -n) + (a 2 - a 2 )(2A -1) + (e + n)( a 2 - A x )] = o. (21) 

The condition (equation (OTI) ) is generally dependent on (j\ and g 2 , but there are two special 
cases. The first is the famous “dark state” | ijj) = -^=(|iV — 1, e, g) — \N — 1, g, e), with the 
condition g\ = g 2 and Ai = A 2 . The spin singlet is decoupled from the photon field, so the 
eigenenergy and eigenstate are coupling-independent. The second case is partly like “dark 
state”—the eigenenergy is also coupling independent, but the eigenstate is not. For gi = g 2 
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and A 1 + A 2 = 1, equation (12TI) reduces to 

(E - N + 1 f[{E — N + l ) 2 + (i — TV)/ - (Ai - A 2 ) 2 ] = 0, (22) 


where g = (j\ +()•>■ For E = N — 1, the condition is ^-independent. Besides, the eigenenergies 
are symmetric about E = N — 1 and there are two degenerate eigenstates with E = N — 1 
existing in the whole coupling regime 


I i>C=N a ) = 
I $C=N b ) = 


1 / 2 (A X — A 2 ) 


A \ y/N — lg 


| A - 2 , e, e) - | N - 1 , e, g) + |N - 1 , g, e) , 


^ ( (A^~ A 2 ) \ N ~ 2 ’ 6 ’ + _ |A- l,^,e) 

Viv^(A 2 -A 1 ) / 


(23) 


(24) 


where -j and 4 are the normalizing constants. For N = 2, i/y;= 2 „) is just \ip) e (equation (fTTl)). 
which is the eigenstate of the two-qubit quantum Rabi model. 

For (7 = 1, the subspace is formed by {|0, e, g), |0, g , e), |1, . 9 , 9 )}, and the eigenvalues 
satisfy 


-E[(Ai — A 2 ) 2 + E(1 — Ai — A 2 ) — 77" + g\ + g 2 \ 

+ (Ai + A 2 — l)(Ai — A 2 )“ + (g\ — g 2 ){Ai — A 2 ) = 0 . (25) 


For gi = g 2 and Ai + A 2 = 1, equation (1251) reduces to 

E[E 2 - X -g 2 - (Ai - A 2 ) 2 ] = 0. (26) 

So there is an eigenstate existing in the whole coupling regime with constant eigenenergy 
E = 0 

\^c=o) = jj ^ Al g A2 ^ |l ,g,g) ~ |0, e, g) + |0, g, e)^ . (27) 

For (7 = 0, the eigenstate is 10, g , g), with constant eigenenergy E = —1. 

To conclude, for identical-coupling gi = g 2 and quasi-resonant condition Ai + A 2 = 
l = c 0 , the spectrum of the two-qubit Jaynes-Cummings Hamiltonian II L]r is very regular and 
interesting: there are horizontal lines at E = N (N = —1, 0,1,...), and the energy curve 
with the same C = N are symmetric about the line E = N — 1. For C = 0,1, there is one 
kind of eigenstates existing in the whole coupling regime with constant eigenenergy, while 
for other cases, there are two such kinds of degenerate eigenstaes, one of which for C = 2 is 
also the eigenstate of the two-qubit Rabi model. With constant eigenenergy, these eigenstates 
are partly like “dark state”, but they are coupling dependent. Choosing A x = 0.7, A 2 = 0.3 
and gi = g 2 = gj 2, the spectra of the two-qubit Jaynes-Cummings model and Rabi model are 
compared in figure [2 
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Figure 2. (a) The spectrum of two-qubit Jaynes-Cummings model with Ai = 0.7, A 2 = 
0.3, w = 1, <71 = < 72 > 0 < <7 = <71 + <72 < 1. (b) The numerical spectrum of two-qubit 
quantum Rabi model with the same parameters. E + and E- are solutions with even and odd 
parity respectively. 


3. Solvability of the two-qubit quantum Rabi model using Bogoliubov operators 

First for convenient, we make unitary transformations S± = -A(a\ x + a iz ) and S '2 = 
(a 2x + a 2z ) to the two-qubit Rabi Hamiltonian (equation O) to obtain (oj is set to 1) 

H' tq = a ] a + gia lz (a + a 1 ) + g 2 a 2z (a + a ] ) + Aiaiz + A 2 a 2x . (28) 

H' tq has a conserved parity with the 7L 2 transformation R = T 0 o\ x 0 a 2x , where T = 
exp(i7ra t a), giving us a way to diagonalize the Hamiltonian in the basis of {|e)i, \g)i}, which 
is the eigenvector of <j\ z . Applying the Fulton-Gouterman transformation HI611401 . 

uJ 1 1 

y T 0 o 2x -T 0 cr 2x 

we obtain 

U '< U = ( 7 H. ) ■ (30) 

where 

H± = a)a + gi(a + a f ) + g 2 (a + a^)a- 2z + A 2 a 2x ± Ai Ta 2x , (31) 

acting on the subspace of R with eigenvalues ±1. First we consider H + . For H_, we just 
need to substitute —for A 1 . In the basis of (|e) 2 0 \4>i), \g) 2 0 \4>2)}, where |0i) and |0 2 ) 
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are photon field states, 77+ is expanded as 

<2+G T C/(a T CJ+) A 2 “1“ AxT 
A 2 T A \T o)a -|- g' (o -|- o+) 


(32) 


To remove the linear terms of a + and a, we use the following Bogoliubov operators 

A = a + g, 77 = cl T g . (33) 

Firstly we use the Bogoliubov operator A. The time-independent schordinger equation reads 

(A+ A — g 2 — E) |0i) + A 2 I 02 ) + Ai|04) = 0, (34) 

[A^A + (</ — g)(A + A+) + S’" — 2gg' — 77]|0 2 ) + ^ 2 |0i) + Ax|</> 3 ) = 0, (35) 

where |</> 3 ) = T|0x), © 4 ) = T\o 2 ). To apply the reflection symmetry, we make the 
transformation T to (l34l) and (l35l) to obtain 

(A+ A — 2 g(A + A+) + 3 g 2 — 77)|</> 3 ) + A 2 I 04 ) + Ai|02) = 0, (36) 

[A^A - (g' + g)(A + A+) + g 2 + 2gg' - E}\<f> 4 ) + A 2 |0 3 ) + Ax|0x) = 0.(37) 

We expand the photon field states \(f>j), j — 1, ..., 4 in terms of the normalized orthogonal 
extended coherent state Ifi2l 

~-9 2 /2-ga t 


K g) = 


^/nJ. 


<a' + g) n , 


which is the eigenstate of A^A, and obtain 

OO 

1 4>j) = e 92/2 ^2Vr^.a jt7l \n,g), j = 


(38) 


(39) 


n=0 


Substituting (l39l) into (l34l)-(l37T). and left multiply (m, g |, we obtain the recurrence relations 


for CLj m 


{E 171 T g )Ox,m — Ax 04,m T A 2 U 2 ,m; 


(40) 


(0' - g)(m + 1)^2, m +1 = (E-m + 2gcf - g 2 )a 2 , m ~ (.</ ~ g)a 2 , m -i 

— Ax a 3 , m — A2ai, m , (41) 

2g{m + l)a 3)m+ i = (m — E + 3 g 2 )a 3jm — 2ga 3>m _x + AxG^m + A 2 a4 jm ,(42) 
(0 + 0')(m + 1)^4,m+i = (m — 77 + 25-5-' + g 2 )a^ m - (g + s0 a 4,m-i 

+ Aiai jm + A2 a 3 , m - (43) 

It is seen the coefficients a^ m depend on three initial conditions, which can be chosen as 

(a+O, ®2,0; 0-3,0}- 

Then we consider the Bogoliubov operator B = a + g'. Now 77+ is given as 

77+77 + (.9 - g'){B + 77+) + («/) 2 - 2.gr/ A 2 + A {T 
A 2 +A {T B^B-g' 


( 44 ) 
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Applying transformation T to the time-independent schodinger equation, we obtain four 


equations similar to (l34l)-(l37l) 

(B^B + (g — g')(B + B^) + ( g ') 2 — 2 gg' — E)\ipi) + A 2 |</? 2 ) + A^^) = 0, (45) 

(B'B - G g') 2 - E)\<p 2 ) + A 2 |<^> + Ail^a) = 0, (46) 

(B'B - (g' + g)(B + B t) + {g'f + 2 gg' - E)\<p 3 ) + A 2 |^ 4 > + A^} = 0, (47) 

(B'B - 2 g\B + B') + 3(g') 2 - E)\<p A ) + A 2 |^ 3 ) + A^) = 0 . (48) 


Expanding the photon field states as | cpj) = e ^) 2 / 2 Yl^=o Vn\bj >n \n, g'), j = 1, • ■ •, 4, where 
the normalized extended coherent state | n,g') is the eigenstate of B, and left multiplying 
(m, g'\. we obtain the recurrence relations for b jm 

(g - 9 ')( m + l)6i,m+i = {E-m + 2 gg' - ( g') 2 )bi,m 

+ (</ — g)bi, m -i — Ai& 4;m — A 2 6 2>m , (49) 

(E — m + (g') 2 )b 2tTn = A\b 3trn + A 2 6 i im , (50) 

(g + g')(m + 1)63, m +1 = (m - E + 2 gg' + ( g') 2 )h, m 

— (g + g')bz,m -1 + A 462,771 + A 2 6 4im , (51) 

2tf'(m + l)6 4 ,m+i = (m - E + 3(g') 2 )fr4,m 

— ‘Zg'b^m-i + A 4 fe 1)m + A 2 b 3tm . (52) 

There are three initial conditions, which can be chosen as { 61 , 0 , b 2< 0 , & 4 ,o}- To utilize the 
reflection symmetry |0 4 ) = T|0 3 ), \(P 2 j = T\o 4 ), finally, we expand the photon states in 
terms of the photon number states as | ipj ) = y/n\cj tTl \n), and obtain the recurrence relations 


for c hm 

{jn T l)^ci jm+ i {E m)c 4)TO gc\ m —\ A 2 c 2 m A 4 c 4 m , (53) 

(m + l)g'c 2jm+1 = (E- m)c 2tm - g'c 2 , m -i - A 2 c 4 , m - Aic 3 , m , (54) 

{m + l)^c 3im+ i = (m - E)c 3:Tn - gc 3 , m _i + A 2 c 4)m + A 4 c 2 (55) 

(m T 1)^ c 4 r7l _)_ 4 (ttt. E)c^ m (7 c 4>m —1 T A 2 c 3 m A \C\ m . (56) 


Considering l^i) = T |^ 3 ), |^ 2 ) = T|^ 4 ), we obtain c hm = (-1 ) m c 3 , m , c 2>m = (-l) m c 4>m , 

so there are only two initial conditions, which can be chosen as ci , 0 and c 2)0 . 

States \4>j), | ipj) and Ty) in different representations should be only different by a 
constant (here can be chosen as 1) if they are nondegenerate eigenstates with eigenvalue E, 
so we obtain 8 equations 

\4>i) = M, (57) 

I <Pj) = iV’j)- (58) 

For practical calculation, we left multiply (0|e^a, where d is chosen arbitrarily, then (1571) 
and (1581) are mapped to 

OO 

< 0 |e fta | <j>j) = Y a j,m e xp( gfii) (/3i + g) m 

m= 0 

00 

= (0|e*“|V’j> = &i,„exp(-s'A)(A + sT, 

m=0 


( 59 ) 
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(o|e^ 2 °| <pj) = Y hm ex p(-g'M^ + g') m 

m =0 

oo 

= (0\e^ j ) = Y c j,m(3?. (60) 

m=0 

Now we are still dealing with power series with infinite terms, so to obtain clear reliable 
result, we must make all the power series convergent. According to the recurrence relations 
for a jrri (equations (I40l)-(l43l)). b jm (equations (l49l)-(l52l)f and c J/rn (equations (I53l)-(l56l)). we 
find the radii of convergence of corresponding power series are \g — g'\, rnirijr; — g', 2 g'} and 
g' respectively. So, for different g and g', we can always choose proper j3\ and d 2 to obtain 
convergent power series m, so that finite terms can give reliable results and by choosing 
proper cutoff, and we can obtain the results with arbitrary accuracy. That is the advantage of 
choosing these three different representations. Because of the linearity of recurrence relations, 
we can denote 


3 

0j(A) = <0|e /3l “| (pj) = Y a k,o0j{f3i), 

k=l 

(61) 


(62) 

k= 1,2,4 


■ft(ft) = {0|eA“|^> = V 6 t , ov >‘(/ h), 

k= 1,2,4 

(63) 

2 

fW 2 ) = (0|e^ 2 “| ^j) 

(64) 

k= 1 



where for example, ^(f3 \) is obtained by setting b k 0 equal to 1 and other initial conditions 
equal to 0 in equations (l49l)-(l52l). like in ||23ll . Now we have eight initial conditions for eight 
equations 

<Aj(/ 3 l) = (65) 

Vjifo) = ( 66 ) 

which can be denoted as 

M jk e k = 0, (67) 

with e = {6i,o, b 2 , o, & 4 ,o, ai,o, « 2 ,o, a 3 , 0 ) ci,o, C 2 ,o} T - The determinant of M, which is just the 
function of energy E must equal to 0, so we obtain 

G + (E) = det(M + ) = 0, (68) 

which can be used to determine the eigenenergy E. Equation (l68l) takes similar form as 
equation (14) in lfT9ll . but are obtained in a simper and more physical way. Choosing Ai = 0.4, 
A 2 = 0.3, u = 1, gi = 3 g 2 , to have convergent power series in equations (l59l) and (l60l) . we can 
choose j3\ = —Sg 2 and (3 2 = —g 2 , then the results can be obtained with arbitrary accuracy. 
The spectrum is shown in figure [3} It is seen there are no level crossings within the same 
parity subspace, so we can label each eigenstates with two quantum numbers—energy level 
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and parity, but the total degrees of freedom are three, so according to the quantum integrability 
criterion proposed by Braak the model is non-integrable, consistent with what the narrow 
avoided crossings in the same parity subspace indicate [l4H and the result in [fT9l . 



g 


Figure 3. The spectrum of two-qubit quantum Rabi model with Ai = 0.4, A 2 = 0.3, w = 
1, g\ = 3 r/ 2 , 0 < <7 = <j\ + <]2 < E n+ and E n _ are numerical solutions with even and 
odd parity respectively, while E+ and E- are analytical solutions with even and odd parity 
respectively. 


4. Conclusions 

We have clarified the algebraic structure behind the possibility of the quasi-exact solutions 
with finite photon numbers found in lfT9l . By analyzing the Hamiltonian structure in the 
photon number space, we find that the permutation symmetry of the qubit-photon coupling 
terms for the two qubits brings about closed subspace, and hence quasi-exact solutions for 
certain parameters. The novel coupling-dependent eigenstates existing in the whole coupling 
regime with constant eigenenergy E equal to single photon energy huj correspond to quasi- 
exact solutions with at most 1 photon, with the condition for the qubits energy splittings 
± A 2 = fku or A 2 — Ai = hu. We have demonstrated this directly from the Hamiltonian 
structure. These special eigenstates are partly like “dark states”, but are coupling-dependent, 
which may have some potential application. Furthermore, based on our study on the two-qubit 
quantum Rabi model, we conjecture such “dark states”-like eigenstates commonly exist in 
similar models with permutation symmetry of the qubit-photon coupling terms. For example, 
for the homogenous coupled two-qubit Jaynes-Cummings model, there are many such kinds 
of eigenstates with constant energy E = Nftw ( N = — 1, 0,1,...) in the whole coupling 
regime, with the condition Ai + A 2 = hu. One of these special states is also the eigenstate 
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of the two-qubit quantum Rabi model. Since the Jaynes-Cummings model is simper than the 
Rabi model, we may find the application of these special eigenstates easier. 

Besides, using Bogoliubov operators, we have analytically retrieved the solution of 
the two-qubit quantum Rabi model. We find three different representations to expand the 
Hamiltonian, and the solutions can be determined by convergent power series. In this way, 
the eigenproblem of the infinite dimensional Hamiltonian can reduces to finite dimensional in 
practical calculation reasonably, and the results can reach arbitrary accuracy. Without using 
Bargmann space, this method is more physical and concise. 
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